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Abstract. This survey paper examines the work of J. von Neumann and 
M.H. Stone as it relates to the abstract theory of wavelets. In particular, we 
discuss the direct integral theory of von Neumann and how it can be applied to 
representations of certain discrete groups to study the existence of normalized 
tight frames in the setting of Gabor systems and wavelets, via the use of group 
representations and von Neumann algebras. Then the extension of Stone's 
theorem due to M. Naimark, W. Ambrose and R. Godement is reviewed, and 
its relationship to the multiresolution analyses of S. Mallat and Y. Meyer and 
the generalized multiresolution analyses of L. Baggett, H. Medina, and K. 
Merrill. Finally, the paper ends by discussing some recent work due to the 
author, Baggett, P. Jorgensen and Merrill, and its relationship to operator 
theory. 



1. Introduction 

This survey paper is an expanded version of a talk given at the AMS Special 
Session on Operator Algebras, Quantization, and Noncommutative Geometry: A 
Centennial Celebration in Honor of J. V. Neumann and M.H. Stone, which was 
organized at the AMS Annual Meeting in Baltimore by Robert Doran and Richard 
Kadison, and took place January 16 - 17, 2003. 

In the mid 1980's, what is now known as wavelet theory was developed into 
a coherent framework, although many of the techniques in this theory had been 
developed as early as the 1960's. This theory involved the analysis of functions 
/ in L 2 (M. n ) by studying the frequency content locally in time. The key theme 
behind this approach is to analyze a function by studying dilates and translates of 
an appropriate family of "wavelets" and/or "Gabor functions", which in turn gives 
either a orthonormal or frame expansion for the function being studied in terms of 
dilates and translates of a finite number of functions. This theory has led to many 
important applications in physics, engineering, and signal processing. In turn, new 
tools in mathematics, in particular a different approach to Fourier analysis, have 
been obtained by analyzing the work of physicists and engineering in subdivision 
schemes and subband coding theory. 
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It is the point of view of some wavelet theorists that to use the Fourier trans- 
form, or more generally, duality for abelian groups, in the study of wavelets, is 
somewhat circular, because one of the points of wavelet theory is to offer an alter- 
native to ordinary harmonic analysis. On the other hand, at certain stages in the 
development of wavelet theory, Fourier analysis has proved critical. This theme is 
illustrated very clearly by O. Bratteli and P. Jorgensen in their recent book jB.Tj . It 
is the intention of this paper to discuss how certain key aspects of operator theory 
due to John von Neumann and Marshall Stone arise very naturally in the study of 
wavelets. Indeed, because wavelet bases and frames and Gabor bases and frames 
come from dilation, translation, and modulation operators acting on a finite set of 
functions, it is natural to expect that the spectral analysis of these operators is im- 
portant in the study of wavelets. Moreover, dilation, translation, and modulation 
operators on L 2 (W l ) are all unitary operators, and various combinations of these 
operators give rise to countable subgroups of the unitary group of a Hilbert space. 

The approach taken in this article, then, is to think of the key operators in 
wavelet theory as arising from unitary group representations of countable discrete 
groups, and then to apply the tools of direct integral theory and von Neumann 
algebras to their study. Direct integral theory in its simplest form is the analysis 
of a single unitary operator, and thus reduces to the spectral multiplicity theory 
of E. Hcllinger and H. Hahn, worked out in 1912-1913. A generalization of this 
theory led naturally to Stone's theorem concerning the decomposition of a unitary 
representation of the real line on a Hilbert space by means of a projection- valued 
measure on K jS] , which in turn gave a decomposition of the original representation 
as a direct integral over K. of one-dimensional (hence irreducible) representations. 
A further generalization of Stone's theorem, due to both Stone and von Neumann, 
led very naturally to a classification of the irreducible representations of the real 

(1 x z \ 
1 y : i,i/,zel}. 
° 1 / 

After a review of direct integral theory for unitary representations, it is our aim 
to survey how work of Stone and von Neumann, specifically in group representations 
and direct integral theory, has been used in the study of the theory of wavelets to 
date. We begin with a discretized version of Stone and von Neumann's approach 
to the study of representations of the Heisenberg group, that is, with a study 
of unitary representations of the integer lattice of this Heisenberg group, H = 
I 1 p r \ 

{ 1 q : p,q,r <E Z} that are closely tied up with the theory of Gabor 
V 1 J 

frames; the representation involved generate translation and modulation operators. 
We discuss a proof due to L. Baggett [15] of a earlier result on such frames due 
to I. Daubechies |D1| and M. Ricffel Ri which uses direct integral theory. We 
then move on to a discussion of wavelets and their construction from the wavelet 
sets of X. Dai and D. Larson |DL| . and discuss a result of L.-H. Lim, J. Packer, 
and K. Taylor jLPTj . In this case the group involved will be a discrete semi-direct 
product group coming from dilation and translation operators on L 2 (R"). Finally 
we move on to multi-resolution analyses and some recent work of L. Baggett, H. 
Medina, and K. Merrill on generalized multircsolution analyses BMM . In these 
results, the ordinary spectral theory applied to unitary representations of Z™ plays 
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a key role. Recent progress in the this area due to Baggett, P. Jorgensen, Merrill 
and Packer will be outlined B JMP . 

Naturally in a short paper of this kind much material must be omitted. In 
particular, we make no attempt to cover wavelet packets, or any of the exciting 
work relating group representations of continuous groups to wavelets. We thank 
Professors Lawrence Baggett, Palle Jorgensen, Wayne Lawton, Kathy Merrill, Ar- 
lan Ramsay, Marc Rieffel, Wai-Shing Tang, and Keith Taylor for many helpful 
conversations on topics related to this paper. 

2. Direct Integrals, Direct Integral Decomposition and Induced 

Representations 

Professor Richard Kadison kindly informed us that the paper of von Neumann 
"On rings of operators. Reduction theory" |vN2| . which concerned direct integral 
theory and its relationship to the decomposition of von Neumann algebras, had 
been written in the late 1930's but sat in von Neumann's drawer at the Institute 
of Advanced Study in Princeton until F.J. Mautner asked him about the material 
in this paper in the late 1940's 1 . Indeed, von Neumann had published results on 
the theory of direct integrals of Hilbert spaces as early as 1938 in jvNl| . and 
he mentions in the preface |vN2j that the bulk of this "paper in the drawer" 
was written during the time period 1937-38. We briefly review the theory in the 
following paragraphs. Although our notation is not the same as von Neumann's 
original notation, the underlying spirit is the same. 

2.1. A review of direct integral theory. We review the theory of direct 
integrals of Hilbert spaces and direct integrals of representations. Here we depend 
heavily on the expositions of J. Dixmier jDxj . R. Fabec |Fa| . G. Mackey |Mac| and 
R. Kadison and J. Ringrose |KR| . 

Let (X, B, v) be a standard Borel space, that is, B is the Borcl structure of a 
Polish topology on a topological space X. 

Definition 2.1. Let X be a Borel space. We say that X is a Hilbert bundle 

over X if there is a surjective Borel mapping p : X — > X such that the following 
conditions are satisfied: 

(i) For each x £ X, the space p~ 1 (x) — Ti x has the structure of a Hilbert 
space. 

(ii) There exists a sequence {fi}^i of Borel functions fi : X — > X with 
P ° fi( x ) — iVi £ X such that the closed linear span of {fi(x)}°l 1 in 
TL X is equal to TL X for all x e X. 

(hi) For all i,j E N, the function x i-^< fi(x), fj(x) >n m is Borel. 
(iv) The function g : X — > X is Borel if and only if p o g is Borel and the 
functions 

x i-» < g(x), fi(p(g(x))) > H p ( g ( x )) 
are Borel for all i. 

The Borel cross-sections to the Hilbert bundle are those Borel functions 
/ : X — > X which satisfy p o f(x) = x. 



Richard Kadison, anecdote given on January 16, 2003, in Baltimore. 
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Definition 2.2. A Borel cross-section / : X — > X is said to be square- 
integrable if J x \\f(x)\\ 2 di' < oo. The set of square-integrable cross sections is an 
inner product space, with the inner product (/, g) given by the obvious formula 

if,g) = < f{x),g(x) > dv. 

Jx 

Upon identifying square-integrable cross sections /, /' satisfying f„ \ \f(x)—f'(x) \ \ 2 dv = 
0, we obtain the direct integral Hilbert space j x \H x ]dv. 

Example 2.3. Let Tt be a fixed separable Hilbert space, and suppose that 
X = X x H. Define p : X — > X to be projection in the first variable. Then 
7i x = 7~C, Vx € H. Fixing a basis {e^} C TL and defining fi(x) = e^Vx G X, we 
see that X is a Hilbert bundle. The vector space of all Borel functions / : X — » 7Y 
such that < f(x),ei > is measurable for every ej gives the Borel cross-sections to 
the bundle. When we consider the direct integral Hilbert space J x \H x ]dv in this 
case, we obtain exactly the family of square-integrable functions from X into 7i, 
where two functions are identified if they are equal except on a set of j/- measure 0. 
The associated Hilbert bundle is called a constant bundle, and the direct integral 
is frequently denoted by L 2 (X,TC) = L 2 (X) ® TL. Specializing further to the case 
H = C, we obtain exactly L 2 (X). 

We note that by condition (ii) of Definition 12.11 each Hilbert space Ti x is sep- 
arable. Define a functions dim : X — > N U {oo} by dim(a;) = dim(7i a; ). For 
each j G N U {oo} let Xj = dim _1 (j). Von Neumann proved that for any Hilbert 
bundle defined as above, and for each j G N U {oo}, Xj is measurable. Hence the 
space X can be partitioned as the disjoint union of measurable subsets 

X = U Xi, 

j£NU{oo} 

and it is known that the direct integral Hilbert space \H,x\dv\x- is isomorphic 
to the direct integral of the constant bundle Ti-(j), where 7i(j) represents a fixed 
separable Hilbert space of dimension j. 

Given a direct integral Hilbert space j x \H x ]dv, we recall how one forms direct 
integrals of bounded operators to act on the Hilbert space. Suppose we are given 
a Hilbert bundle X, and suppose that x <— ► T{x) is a function that assigns to each 
x € X an element of B(Ti x ) 1 the algebra of all bounded operators on Tt x . 

Definition 2.4. Let (J,p,X) be a Hilbert bundle. We say that the mapping 
x i ► T(x) is a measurable cross-section of bounded operators if for every Borel cross- 
section f : X — > X, the cross section x i— > T(x)f(x) is also in X. If x i— * T(x) is 
a measurable cross-section of operators, we say it is a bounded measurable cross- 
section if the function x i— > ||T(x)|| is in L°°(X). In this case, x >—>■ T(x) determines 
an element of B(f®\H x \dv), denoted by 

T = T [T{x)]du, 
Jx 

where we have T(/)(x) = T(x)(/(x)) for any square-integrable cross sections 
/ € [H x ]dv. Conversely, we say a bounded operator T € B(f x [H x ]dv) is 
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decomposable if there is a bounded measurable cross-section x i— > T(x) such that 

T = [T{x))dv. 
J x 

Now we briefly discuss the direct integral decomposition of unitary group 
representions. Suppose that G a locally compact second countable group. Let 
U : G — ► IA{TL) be a continuous unitary group representation of G, that is, a 
homomorphism of the group G into the unitary group U(TL), such that Vx, y € 7i, 
the map 

g i-> < U g (x),y > 

is continuous from G to C. 

Definition 2.5. We say that Z7 is decomposable as a direct integral over the 
Borel space (X, v) if the following conditions hold: 

(i) The Hilbert space H is isomorphic to a direct integral j^\H x ]dv. 

(ii) For each x € X there is a unitary group representation U{x) : G — > H x 
such that 

U g = f [U(x) g ]du, V.9 e G, 

with respect to the identification of Tt with J x \H x \dv made in (i). 

The theory of von Neumann algebras and in particular, the paper |vN2) . was 
fundamental in showing that every unitary decomposition of G could be decom- 
posed as a direct integral of irreducible representations of G. In order to briefly 
describe this procedure, we first recall the definition of a von Neumann algebra 
M of B(H), due to F. Murray and J. von Neumann jMvNll IMvN2| . As other 
articles in this survey article will go into these definitions in more detail, here we 
deliberately make our remarks sketchy and refer the reader to |Dxj . for example, 
for full details. 

Definition 2.6. A von Neumann algebra on the Hilbert space Ti. is a sub- 
algebra M. of B{TL) that is closed under the adjoint operation, contains the identity 
operator Id-n, and is closed in the strong operator topology. The set of all operators 
in B(TL) that commute with every element of M. is called the commutant of M. 
and is denoted by M! . 

Murray and von Neumann proved that if AW is a von Neumann algebra, then 
its commutant M! is a von Neumann algebra also, and 

[M']' = M. 

Indeed, if N is any subalgebra of £>(7i) which is closed under the adjoint operation, 
they proved that the smallest von Neumann algebra in B(7i) containing M is . 
For this reason the double commutant JV" is called the von Neumann algebra 
generated by Af. A von Neumann algebra M. is said be a factor if 

M n M' = {Aid : A e C}. 

This means exactly that the center of M. consists of scalar multiples of the identity 
operator. Von Neumann and Murray showed that one can define a dimension 
function, which is unique up to a normalizing scalar, on projections of a factor M., 
which extends to a trace function on the factor M.. A factor M. is said to be of 
type Hi if this dimension function takes on exactly the values [0, 1]. 
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Of course not all von Neumann algebras are factors, and there is a vast literature 
describing the structures and decomposition of von Neumann algebras. Here we 
only touch on that part of von Neumann algebra theory relevant to our discussion. 
We first use von Neumann algebras to describe how to decompose unitary repre- 
sentations of groups into direct integrals. This procedure is due to von Neumann 
|vN2| . F. Mautner [Mauj and R. Godement |Goj . and goes as follows: 

(i) Form Qu = {U g : g £ G}" , the von Neumann subalgebra of B(TL) 
generated by {U g : g £ G}; 

(ii) Let M. be a maximal abelian subalgebra of the commutant of Qu, and 
find an isomorphism between M. and L°° (X, v) for some measure space 

(iii) Decompose Ti as J x \H x ]dv and U as J x [U(x) g ]dv. 

However, since the maximal abelian subalgebras in Qu' are not unique up to inner 
conjugacy, this direct integral decomposition need of U not be unique, and in the 
case where G is countable and discrete and [G : [G, G]] is infinite, there always exist 
representations having non-uniqueness in their direct integral decompositions; see 
|Macj . Chapter 3, Section 5 for an example. 

2.2. A review of induced representations. Induced representations will 
play an important role in the discussion that follows, so we review their definition 
at this point. We follow the description of induced representations given in Chapter 

x of IE]. 

Definition 2.7. Let G be a locally compact second countable group with 
closed subgroup N, and suppose we are given a Borel cross section s : N\G — > G 
with s([eo]) = eg, where here e<3 represents the identity element of the group G. 
Use s to define a one-cocycle a for the right action of G on the coset space N\G 
as follows: a : N\G x G — > N is defined by 

a(x,g) = s{x)g[s{x ■ g)]- 1 , (x,g)£ N\G x G. 

Let p be a unitary representation of the group N on the Hilbert space TL. Then 
the induced representation Ind,Q (p) has as its representation space L 2 (N\G,Tt). 
where the measure v on N\G is quasi-invariant under translation by G. If this 
measure is invariant, then the formula for Ind ( fr(p) is given by 

Ind%(p)(g)f](x) = p(a(x,g))[f(xg)}, 

for g £ G, x £ N\G, f £ L 2 (N\G,H). 

In the discussion that follows, all of the groups we study will be discrete, and 
the translation-invariant measure on the coset space N\G will just be the counting 
measure. 

Example 2.8. Let H denote the integral Heisenberg group; that is, H consists 
of all matrices of the form 

/ 1 p r \ 
{ 1 q : p,q, r £ Z}. 
V 1 J 

It is easier notationally to parametrize H by the set {(p, g, r) : p, q, r £ Z}, with 
obvious operations coming from the matrix multiplications. Let N be the abelian 
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subgroup of H denned by N — {(0, q, r) ; g,r£ Z}. Fix t, u G [0, 1) and define a 
one-dimensional representation or character x( u > of on C by 

x(«,*Wr) = e 2 ^V™ rt , (0,g,r) G AT. 

The coset space Hq\H can be identified with Z, and the cross-section s : N\H = 
Z — > is given by 

s(fc) = (0,0, fc), fc G AT\H. 
The corresponding one-cocycle for the action of _ff on N\H is given by 

(p, g, r)) = (p,q,0). 
Hence Ind%(x(u,t)) is defined on ^ 2 (Z,C) ^ ^ 2 (Z) by 

Ind»( X (u,t))(p,q,r)f(k) = e 2 ™« u e Mrt f (k + r) , fel 2 (Z), (p, q, r) G ff. 

We call the representations Ind^(x(u,t)) "Stone-von Neumann" representa- 
tions of H, because they can be extended to the three-dimensional real Heisenberg 
group and for t ^ provide irreducible representations of the real Heisenberg group 
of the type first discovered by Stone and von Neumann. We remark that elsewhere 
in the literature these representations are referred to as "Weyl-Heisenberg" rep- 
resentations, since W. Heisenberg first described commutation relations between the 
position and momentum operators, and H. Weyl, using the work of von Neumann 
and Stone, noticed that by exponentiating these operators, one could construct 
unitary operators on a Hilbert space satisfying a key commutation relation that 
will arise later in this paper. 

We note that the operations of taking direct integrals of representations and 
inducing representations commute: That is, if N is a closed subgroup of the locally 
compact second countable group G, and if L is a representation of N which is 
written as a direct integral 

/■© 

L n = [L n (x)]dv, ne N 
Jx 

on the Hilbert space J® \H x \dv, then the induced representation Ind^j{L) is uni- 
tarily equivalent to the direct integral of representations 

/ [Ind%{L x )]dv 
Jx 

defined on the direct integral Hilbert space 

lL 2 (N\G,H x )]dv. 

c 

We refer the reader to jMacj and |Fa| for details. 

3. Important operators and groups for wavelet and frame theory and 
the von Neumann algebras they generate 

We first review the notion of a normalized tight frame in a Hilbert space. 

Definition 3.1. A sequence {/„ : n G N} of elements in a Hilbert space 7i is 
said to be a frame if there are real constants C, D > such that 

oo oo 

(3-1) C]T|(/,/„)| 2 < Il/H 2 < D^KfJn)] 2 

n—1 n—1 
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for every / S TL. The constants C and D are called the frame bounds for the 
frame. If C = D, the frame is called a tight frame, and if C — D — 1, the frame is 
called a normalized tight frame, abbreviated NTF. It is known that {/„ : n e N} 
is an NTF for TL if and only if the following reconstruction formula holds: 

oo 

(3.2) / = /„>/„, V/eW. 

n=l 

Because of this last identity, called the reconstruction formula, normal- 
ized tight frames are sometimes referred to as "Parseval frames" in the literature. 
Clearly, an orthonormal basis for TL is a special case of a normalized tight frame. 

3.1. Key operators in Gabor and wavelet frame theory. It turns out 
that in the construction of Gabor and wavelet frames, there are three key types of 
unitary operators that play important roles. They are dilation operators, transla- 
tion operators, and modulation operators, which are defined on L 2 (M"). 

Definition 3.2. Let A be an n x n integral matrix all of whose eigenvalues 
have modulus greater than one. The dilation operator corresponding to A is defined 
on L 2 (R n ) by 

D A (f)(t) = \det(A)\Kf(At). 

Definition 3.3. Let T be a cocompact lattice in K™, for example, F = Z", and 
let 7 € r. The translation operator corresponding to 7 is defined on L 2 (W l ) by 

T 7 (/)(i) = /(t- 7 ). 

Definition 3.4. Let a <G R™. The modulation operator corresponding to a 
is defined by 

M a (m) = e^-'fit). 
Let T denote the Fourier transform on £ 2 (R), defined by 

nf)(?) = I f{t)e~ M dt. 
Jr 

Then one can verify the following relations: 

T 7 D A = D A T A ^, 

M a T fj = e-^^TpMa, 
TDaT- 1 = D AT -i, 

and 

TT^T~ X = M 1 . 

We now use the operators listed above to define certain kinds of frames in 
L 2 (R n ). 

Definition 3.5. Fix full-rank lattices C and K in K™. By this we mean that 
C = L(Z") and JC = K(Z n ) for two non-singular n x n real matrices L and K. We 
say that a vector / e L 2 (R n ) forms a Gabor type NTF for L 2 (R n ) (respectively, 
a Gabor orthonormal basis for L 2 (K™)) with respect to the lattices £ and /C if the 
set of functions 

{M x T K {f) : X e C, k e JC} 
forms an NTF for for i 2 (R") (respectively, a orthonormal basis for L 2 (R")). 
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Specializing to the case n = 1, we say that a vector / £ L 2 (M.) forms a Gabor 
NTF for L 2 (M.) with respect to the constants a, j3 > if the set of functions 

{M pa T Pq {f) : p,q £ Z n , 1 <im} 

forms an NTF for for L 2 (R). 

Example 3.6. Let n = 1 and a /3 = 1, and take f(t) = X[o,i](t)- An easy 
calculation with Fourier series shows that {M q T p (f) : p, q 6 Z} is an orthonormal 
basis for L 2 (R), hence is a Gabor type NTF. 

Example 3.7. Again take n — 1 and a (3 — 1. Let / be the Gaussian 
/(t) = e~* , then the set of functions {M q T p {f) : p,q £ Z} form a set whose finite 
linear combinations are dense in L 2 (K). Von Neumann knew this was true in the 
1940's. However, in this case {M q T p {f) : p, q £ Z} is not a frame, c.f. |Benj . [Pi] . 

Definition 3.8. A NTF wavelet family for dilation by A is a subset 
{■01, • • • , C L 2 (R") such that 

{D^Oi) : l<KmjeZ,«eZ n } 

is an NTF for L 2 (R"). If the functions {D^T^ipi) : 1 < i < m, j £ Z, v £ 
Z n } form an orthonormal basis for L 2 (W l ), the family {ipi, • • • , V'm} is called an 
orthonormal wavelet family for dilation by A. 

Thus when one has a NTF wavelet family {ipi, • • • , ip m } for dilation by A, for 
every / € L 2 (M n ) the reconstruction formula Eauation l3.2l is satisfied with respect 
to the NTF {D j A T v (ip t ) : 1 < i < m, j £ Z, v E Z n }. 

In applications, the function / being studied in L 2 (R") determines the wavelet 
family {tpi, • • • , ip m } selected, and one would hope that all but finitely many of the 
wavelet coefficients in the reconstruction formula for / would be zero or at least so 
small as to be safely ignored in computations, without loss of essential data. 

Example 3.9. Let n = 1 and A = 2, and define 

f 1 if x£ [0,i], 
i/> H (t) = I -1 if x£ (i,l] 
I otherwise. 

Then {D j 2 T v {i[)h) ■ j,v £ Z} forms an orthonormal basis for L 2 (R), and ipn is 
called the Haar wavelet for dilation by 2. 

3.2. Group representations associated to Gabor and wavelet theory. 

The operators connected to Gabor frames and wavelet frames are very much con- 
nected to certain groups, as we describe now. 

Definition 3.10. Let H denote the integral Heisenberg group defined in Ex- 
amplejH Define U°f q r) (/)(*) = e - 2nira P e 2mqat f{t-pj3). Then U a ^ is a unitary 

representation of H on L 2 (K). A simple change of variable in t shows that U a ^ is 
unitarily equivalent to U al3,1 so depends only on the product a/3 and not on the 
separate parameters a and (3. We write II s for U 5 ' 1 . 

Note that 

f/ 5 (p,0,0) = T p , 

and 

U S {0,q,0) = M qS . 
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It follows that / € L 2 (M) gives an Gabor NTF (respectively, a Gabor orthonormal 
basis) for L 2 (R) if and only if {U s (p, q, 0)(/) : p, q £ Z} forms an NTF (respectively, 
an orthonormal basis) for L 2 (R). This description can be extended full-rank lattices 
in the 2n+l dimensional real Heisenberg group as well, and we leave the calculations 
to the reader. Thus the study of Gabor frames is intimately connected to the study 
of the representation II s of H and to the von Neumann algebra generated by the 
operators {U d (p, q,r) : (p,q,r) £ H}. In turn, these representations of H can be 
viewed as projective (or "multiplier" representations) of Z 2 corresponding to the 
multiplier tog on Z 2 defined by 

w«((pi,gi),(P2,«0) = e-*** 1 *', (pi,qi),(p2,to) e z 2 . 

Such representations have been studied extensively in the literature, c.f. |BKURie| . 

We now discuss discrete groups essential to the study of wavelets. We note that 
for the classical case where n — 1 and A = 2, it was X. Dai and D. Larson and, 
independently, B. Brenken, who first realized the connection between the study of 
representations of this group and wavelet theory, c.f. [DLL IBr) . 

Definition 3.11. Let A be an integral nxn dilation matrix, and let Qa denote 
the discrete A-adic rational subgroup of Q" defined by 

Qa = U je z[A- j (Z n )}. 

The matrix A defines an action of Z on by automorphisms which we denote by 
9 : 

9(h)(0) = A- k ((3), /3 e Qa, keZ. 
We thus can form the semidirect product of Qa and Z, Qa X0 Z, which we call the 
wavelet group (for dilation by A) and denote by G a- 

We parametrize Ga by pairs {(/?, k) : (3 S Qa, k G Z}, where the group 
operation is defined by 

(PiM) ■ (f3 2 ,k 2 ) = (/3i + 0(fci)(/3 2 ),fci + fc 2 ). 
Ga is sometimes called a generalized Baumslag-Solitar group, c.f. |MV| . 

Definition 3.12. Define a group representation W of the wavelet group Ga 
on L 2 {WL n ) by 

WW) (/)(*) = \det(A)\^f(A k (t-(3)), ((3,m)eG A . 
We call this the wavelet representation of Ga- 
One calculates that 

W {m = T P D\. 

From this observation we see that {ip%, • • • , tp m } is a NTF wavelet family (respec- 
tively, orthonormal wavelet family) for dilation by A if and only if 

: P e e Z, 1 < i < m} 

forms a NTF (respectively, an orthonormal basis) for L 2 (W l ). Again it follows that 
the study of wavelets is intimately connected to the study of the representation W 
and the von Neumann algebra generated b the unitary operators {W^ ^ : (f3, k) G 
G A }- 

Much of the remainder of the paper will be devoted to the decomposition of 
the representations II s and W via direct integral theory. 
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4. The use of direct integrals in the study of Gabor frames and 

wavelets 

We now use the representations described above in the study of properties of 
group representations and wavelets. 

4.1. The study of Gabor frames using direct integral theory and von 
Neumann algebras. Fix 8 > 0, and consider the Gabor representation U s 
defined in the previous section. The following decomposition theorem is explicitly 
due to L. Baggett, although various people had been aware of it previously; we 
discuss higher dimensions at the end of this section. 

Theorem 4.1. B Let H and N denote the integer Heisenberg group and its 
doubly generated abelian subgroup N — {(0, q, r) ; q, r G Z} defined in \2.8l Let E 
denote the interval [0, 8) and let v denote Haar measure on E. 

(i) There is a unitary equivalence 

U : L 2 (R) -> / [l 2 {1)]dv 

J E 

and a one parameter family of representations {U s (t) : t G E} of H 
on / 2 (Z) such that U s is unitarily equivalent via IA to the direct integral 
representation 

/V(*)]*' 

J E 

of H on f® [I 2 (Z)]dv. Fort G E, the representation U s (t) is defined by 
U s (t) {p ^ r) g(k) = e-^e-^e^Xgik-p), gtl 2 (Z). 

(ii) For each t G E, let xt.s denote the one- dimensional representation of N 
on C defined ro lff.&l i.e. 

XM ((0,9,r)) - e-^e- 27 ^, (0,g,r) G N. 
Then each representation U s (t) of H is unitarily equivalent to the mono- 
mial representation Ind N (xt.s) of H on l 2 (N\H) induced from the one- 
dimensional representation Xt.s of N on C. 

It follows from the above theorem that the study of the monomial Stone-von 
Neumann representations U s (t) of H are key to understanding Gabor frames. It is 
possible to use this direct integral decomposition to approach the following question. 
Question: Does there exist a single function / G L 2 (M.) such that the family 

{/„,,(«) = e 2 ^ at f(t-pf3): p,q G Z}, 
or equivalently, the family 

{K^o)(f) ■■ p,q e Z}, 

forms a NTF for L 2 (R)? 

Answer of I. Daubechies |Dlj : (1) If S = 1, the answer to the question is yes 
fExample 13.01 gives an orthonormal basis {f p , q } )', (2) if 6 < 1, the answer is yes, 
(3) if 6 > 1, the answer is no. 

The first proof of (3) for S > 1 and rational is due to Daubechies |D1| . The 
proof of (3) for 8 > 1 and irrational is due to Rieffel who used von Neumann algebra 
"coupling constant" techniques described in |Ri| . He related this work to his study 
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of the irrational rotation algebra in ,Rie , as certain representations of H give rise 
to this well-known C*-algebra. We sketch Rieffel's proof in later paragraphs, but 
we first describe how this problem has been studied by using direct integral theory. 
In 1990 L. Baggett gave a proof of the same results involving the explicit direct 
integral decomposition of the Stone- von Neumann representation mentioned above. 

Example 4.2. We briefly discuss the above result from the direct integral point 
of view: for simplicity, we consider L 2 (R) within the framework of Baggett's direct 
integral decomposition 

/•© 

L 2 (R) = / [l 2 {Z)]dv, 

J [0,5) 

and consider for what values of S the function 

fs = X[o,S) = / [&o]dv 

J [0,5) 

can generate a tight frame under the orbit of the operators 

{F U\t) M) dv}. 

J [0,6) 

(Here 5q represents the Dirac delta function at defined on Z.) 
We note that 

{UL q ,o)(fe):p,qeZ} = {F U s {t) {p , q , 0) f s {k)dv:p,q,k£l} 

J [0,5) 

= { / e- 2mqt e 2wlqSk f s {k-p) :p,q,keZ} = { e - 2 " qt e 2nlqSk S p (k) :p,q,ke Z}. 

J [0,6) J [0,6) 

When we consider the direct integral as a tensor product L 2 [0, S) <8> ? 2 (Z), this last 
set becomes 

{e 2ntqSk e~ 2jr%qt ®S p : P ,qeZ}. 
Let 6 < 1. From our above discussion, one sees that {U/ pq0 \(fs) '■ p,q £ 2} will 
form a frame for L 2 (M.) if and only if the family of functions {e 2mqt : q £ Z} forms 
a frame for L 2 [0, 5). Let us consider the normalized family {-j=e 2mqt : q £ Z} for 
5 < 1. We take / £ L 2 [0, S) and compute 



EK/.4 



2iriqt\ rftr r j2 



■■e )L*[0,6)\ 



qei, 

(since / is supported on [0, 6)) 



(by Parseval's formula) 



J2\ f f(t)^e^Hu(tf 

J [0,6) VO 



[0,5) 

f(t)e- 27rlqt dv{t)\ 2 

[0,1) 



tII/III 2 [o,i) 



jll/lli 2 [0,5) 
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(since L 2 [0, 8) C L 2 [0, 1)). Thus if < 8 < 1 we obtain a tight frame for L 2 [0, 8) <g) 
£ 2 (Z) with frame bound 8. Indeed if 8 = 1 we get an orthonormal basis via this 
method, as we have already seen in Example 13. 61 

On the other hand, if <5 > 1 the set {e 27 ™ 9 * : q € Z} does not even span L 2 [0, S) : 
so one cannot hope to get a frame from {U^ p -\(fs) ■ Pi Q e ^} m this case. 

We now discuss this result from the point of view of M. Rieffel, who gave the 
first proof that if a ■ /3 was irrational and greater than 1, it was impossible for a 
vector / 6 L 2 (R) to form a Gabor NTF in L 2 (M.) with respect to the constants 
a, 0. |Rij . As before, set 8 — a ■ (3. A special case of a result of M. Takesaki 
Ta shows that is possible to form the finite transformation group von Neumann 
algebras M(8Z,R/Z) acting on L 2 (R) on the left and M(Z, (<5Z)\R) acting on 
i 2 (M) on the right in such a way that each von Neumann algebra is the other's 
commutant. [We note that it is now standard to write such transformation group 
von Neumann algebras as L°°(R/Z) x 8Z and L°°((<5Z)\R) x Z respectively]. Rieffel 
showed that M(Z, (<5Z)\R) is generated by two operators U and V which satisfy 
the standard commutation relation for the irrational rotation algebra 

UV = e 2mS VU. 

Rieffel's work showed that there was a unitary equivalence W : L 2 (R) — > L 2 (R) 
such that 

Wtf?i,o,o)W = V 

and 

W^fo,i,o)W* = U. 

It then became clear that in order that there exist / G L 2 (M.) such that 
{U? p q o)(/) : P, 9 G Z} forms a NTF for L 2 (R), it was necessary that the von Neu- 
mann algebra M(8Z, R/Z) should contain a cyclic vector in L 2 (R). Now M(8Z, R/Z) 
and M(Z, (<5Z)\R) are both finite von Neumann algebras in the sense that they have 
defined on them canonical center-valued traces, which shall be denoted here by r 
and r', respectively. Indeed, when 8 is irrational, these von Neumann algebras are 
in fact Hi factors, so that we can think of the traces as being scalar multiples of the 
identity. Thus in the case of irrational 8, we can identify the range of these traces 
with complex numbers, and in the case where the traces are defined on projections 
in the respective von Neumann algebras, the traces will take on real values between 
and 1. Fixing some notation, suppose that M and its commutant M' are two finite 
von Neumann algebras acting on a Hilbert space Tt. Then for £ £ Ti. let be that 
projection in the center of M corresponding to the smallest closed subspacc of H 
containing M'(£), and let E'^ be that projection in the center of M' corresponding 
to the smallest closed subspace of 7i containing We remark that if £ is a 

cyclic vector for M, then E'^ = Wm' ■ Recall from |Ri| that the coupling function 
lu is the possibly unbounded operator defined from the center of M to the center 
of M' which satisfies the condition 

r'(E' e ) = rH^)). 

Now suppose that M(5Z,R/Z) = M and M(Z, (<5Z)\R) = M' In Theorem 3.2 and 
Example 4.1 of |Ri| . Rieffel computed the coupling function to corresponding to the 
pair (M,M r ) in this case to be exactly multiplication by the positive number 4. It 
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follows that if / is any vector in L 2 (R), 

t'(E>) = yr(E f ), 

which is strictly less than one if S > 1 (since < r(Ef) < 1). As t'(E'j) < 1, E'j 
will not equal to Id^/ ' , which by our above remarks shows that / cannot be cyclic 
for M. 

In addition to the work of I. Daubechies, M. Rieffel and L. Baggett on the 
Gabor frame problem (sometimes called "Weyl-Heisenberg frames" in the litera- 
ture), Daubechies, H. Landau, and Z. Landau |D3| also pursued different methods 
in dimension greater than 1, as did J. Ramanathan and T. Stegcr RS]. Both of 
the papers cited above involved the use of von Neumann algebras in their proofs. 
A. Ron and Z. Shen |RoS) also worked on this problem using a different approach. 
Very recently, D. Han and Y. Wang |DL| have proved the following generalization 
of results on Gabor frames, which applies to lattices in R ra , for n > 1. 

Theorem 4.3. HW Let n G N, and let C and JC be two full-rank lattices in 
R", called the modulation and translation lattices, respectively. For A G C, k G /C, 
and g G L 2 (R"), let 

9\.A X ) = M\T K g(x), 
where M\ and T K are the standard modulation and translation operators. Then 
there exists g G L 2 (W n ) such that {<7a,« : A G £, k G K.} form a frame for L 2 (R") 
if and only if v(C)v(IC) < 1. 

Here v(C) and v(IC) denote | det L\ and | det K \ respectively, where L and K are 
n x n matrices with C = L(Z n ) and K, = K(Jj n ). We remark that the necessity of 
the volume condition had been proven in |RSj and |RoS) and it is the sufficiency 
that was established in |HW] . The proof of Han and Wang uses geometric lattice 
tiling and packing methods. 

It would be interesting to approach the higher dimensional case by decomposing 
a certain representation of a discrete subgroup of the 2n+ 1-dimcnsional Hcisenbcrg 
Lie group as a direct integral. 

4.2. The study of wavelets using the direct integral of the wavelet 
representation. Recall Ga denotes the discrete semi-direct product group QU x g 
Z, where A is an n x n integral dilation matrix and denotes the discretized 
A-adic subgroup of Q n obtained from A, viewed as a discrete group. For example, 
if n = 1 and A is the constant dilation by 2, Q2 represents the group of dyadic 
rational numbers. 

Recall from Definition 13. 121 that the wavelet representation W of Ga is given 
by the formula 

W {m = T p D k A , (f3,k) G G A . 

We discuss how to decompose this representation as a direct integral of represen- 
tations. The notion of wavelet set due to X. Dai and D. Larson in the classical 
A = 2 case DL|, and for more general dilation matrices to X. Dai, D. Larson and 
D. Specgle DLS is needed. 

We consider dilation and translation in the frequency domain, using the nota- 
tion of Section 3.1. Denote 

TJ2 = TDaT- 1 , % = TT V T~ X , v G Z". 
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Definition 4.4. f |DLL IDLSj ) Let A be an n x n integer dilation matrix on 
R™, with A T denoted by B. We say that a set E C R™ of finite positive measure 
is a wavelet set for dilation by A if , 1 ye is a wavelet for dilation by A in the 

frequency domain, that is, if 

(d^t v {—1=xe) ■ j e z, v e i) 

is an orthonormal basis for L 2 (R"). 

This definition would imply by use of the inverse Fourier transform that 
J-~ 1 ( , 1 ye) is a wavelet for dilation by A. 

Dai and Larson have shown that E C R™ is a wavelet set for dilation by A if 
and only if the following conditions are satisfied: 

(i) E is a measurable subset of R n , with B j (E) H B k (E) = 0, j ^ k £ Z, 

(ii) KR"\U ieZ I? J '(£))=0; 

(hi) E is translation congruent to the set [—\, \) n modulo the lattice Z n . 
Conditions (i) and (ii) are equivalent to E tiling R n under dilation by B, and 
condition (Hi) is equivalent to E tiling the plane under translation by Z™. 

The following theorem of X. Dai, D. Larson, and D. Speegle show that wavelet 
sets always exist. 

Theorem 4.5. |DLS| Let A be an n x n integral dilation matrix. Then there 
exists a wavelet set in R™ for dilation by A. 

Example 4.6. [DL] Let n = 1, A = 2 and set E = [-1, -1/2) U [1/2, 1). It is 
clear that E tiles the plane by dilation by 2 and by translation by Z. Hence E is a 
wavelet set and ^>(x) — xe(%) is a wavelet in the frequency domain for dilation 
by 2. The corresponding wavelet in the time domain is the Shannon wavelet related 
to the sine function. 

We now use wavelet sets and direct integral theory the wavelet representation. 
We first consider a representation equivalent to the wavelet representation. 

Definition 4.7. The representation W :: (QU xZ-> U(L 2 (W 1 )) defined by 

W((3,k) = TW(f5,k)T- x 

is called the wavelet representation of the wavelet group Ga in the fre- 
quency domain. One calculates that 

W(P,k)(f)(x) = \det(A)\-ie 2 ^ x f(B k x), {J3,k) £ Ga, f £ L 2 (R n ). 

Now let £ be a wavelet set for translation by A; by Theorem 14.51 such a set 
will exist. Note that since E tiles R™ by dilation by powers of B, it is possible to 
identify R™ with the Cartesian product space E X Z, and hence we can set up a 
unitary operator carrying L 2 (R") to L 2 (E xZ)~ L 2 (E,l 2 (Z)), where E carries 
ordinary Haar measure. Upon doing this, and studying the wavelet representation 
of Ga in the frequency domain, we obtain the following result: 

Lemma 4.8. I'i.l'A |LPT] ) Let A be a nx n integral dilation matrix and let E 
be a wavelet set for dilation by A. The wavelet representation W of Ga is equivalent 
to the representation W : Ga —> U(L 2 (E x Z)) defined by 

W i0 , k) (f)(x,m) = e- 2mAm ^f(x,m-k), (0,m) £ Q A , f£L 2 (ExZ). 
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The lemma above allows us to decompose the wavelet representation of Ga 
as a direct integral of monomial representations, just as Baggett did earlier in the 
Gabor case. The following result is joint with L.-H. Lim and K. Taylor LPT : 

Theorem 4.9. |LPT| Let A be a n x n integral dilation matrix, and let E be 
a wavelet set for dilation by A in the sense of \4-4\ Then: 

(i) The wavelet representation W of Ga is equivalent to a direct integral 
of representations J® W x dv of G a on the direct integral Hilbert space 

(ii) Each W is a representation of Ga on l 2 (1i) which is induced from a char- 
acter Xx on the subgroup <Qa of G. If x ^ 0, W x is irreducible. 

We briefly discuss the proof of the above theorem. Firstly, the discussion on 
direct integrals in Section 2 shows that the Hilbert space L 2 (E x Z) can be iden- 
tified with L 2 (E,l 2 (Z)) which in turn can be identified with the direct integral 
Hilbert space J E [l 2 (1)]dv. Upon making this identification, it becomes clear that 
the representation W can be decomposed as a direct integral in the desired fashion. 
The fact that each W x is induced from a character on Qa follows from calculations 
done in |LPT| . The irreducibility of the representation W x for x ^ follows from 
the observation that the orbit of Xx under the induced action of Z in the compact 
abelian dual group Qa is dense whenever x =/= 0. 

Before making some more observations, we recall the definition of the (left) 
regular representation of a group G. 

Definition 4.10. Let G be a second countable locally compact group. The 
(left) regular representation of G is the representation A of G on L 2 {G) given by 

\(f)(h) = fig-'h), g,h eG, f G L 2 (G). 
(Here G is given its left-invariant Haar measure). 

We remark that in their original study of wavelet sets |DL| . Dai and Larson 
actually computed the commutant of the wavelet representation of Q2 x Z and 
showed that it was a commutative von Neumann algebra. They also noted that 
the von Neumann algebra generated by the regular representation of Q2 * ^ was 
a /ii-factor, as was its commutant. Hence the wavelet representation and the 
regular representation of Q> 2 * Z cannot be unitarily equivalent. The same argu- 
ment shows that similar statements are true for the regular representation and the 
wavelet representation of <Qa ^gZ,. By our discussion in Section 2, this implies that 
there is up to unitary equivalence, only one direct integral decomposition of the 
wavelet representation W of Ga, whereas using the method outlined in |Macj . one 
can decompose the regular representation of Ga into two distinct direct integral 
decompositions which are not equivalent. 

One can use Lemma FQ1 to prove the following Corollary, which in the classical 
case is due to F. Martin and A. Vallette |MVj . We first recall that two unitary 
representations of a locally compact group G are said to be weakly equivalent 
if the kernels of the associated representations of the group C* -algebra C*{G) are 
equal. 

Corollary 4.11. |LPT| The wavelet representation of Ga is weakly equivalent 
to the regular representation of Ga- 
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We end this section by discussing an example which shows how to construct 
tight frames in the classical case, when n — 1 and A is dilation by 2, by using the 
direct integral approach. This example is analogous to the approach of Baggett in 
the Gabor frame case. 

Example 4.12. Fix a,b > and let E = [-2a, -a) U [6,26) C R. Let 
A = (2); note that dilation of £ by 2 tiles the real line R. We claim that 
ip = T~ x (\j ' \J h(E)xe) is a tight frame wavelet for translation by the integers and 
dilation by 2 if and only if a+ b < 1. If a+ b < 1, the vectors {D^T k if; \m,k eZ)} 
form a tight frame for L 2 (R) . If a + b = 1 ip is an orthonormal wavelet for dilation 
by 2. 

We do not do the calculations in detail, but note the main idea for this calcu- 
lation is very similar to that given in Example 14. 21 for the Gabor frame case. We 
begin with the case a + b = 1; then setting b = 1 — a, it is already known (c.f. jDLj . 
Example 4.5) that [—2a, —a) U [1 — a, 2 — 2a) is a wavelet set for dilation by 2, and 
the statement follows. If a + b < 1, E still tiles R under dilation by 2. Using the 
methods of Lemma 14.81 it is possible to show that the wavelet representation is uni- 
tarily equivalent to a representation over a direct integral Hilbert space J® [I 2 (Z)]aV. 
Using this representation, the functions {D™T k t\} : m,k £ Z} correspond to the 
functions {-^=e^ tkx <g> S m : k £ Z} C L 2 (E) (g> l 2 (Z) ~ J® [Z 2 (Z)]aV. Since 

{5 m : m £ Z} is an orthonormal basis for / 2 (Z), it follows that {D™Tkip; m, k £ Z} 
will form a frame for L 2 (R) if and only if the functions { , 1 e -^ikx ^ g form 

a frame for L 2 (E). If a + b < 1, one computes that {—j==e~ 27rlkx ; k £ Z} form a 

tight frame for L 2 (E) with frame constant v(E) = a + b just as in Example 14.21 

Finally, if a + b = v{E) > 1, the closed span of the set {1/ \/ )i(E)e~ 2mkx | k £ 
Z} in L 2 (E) will consist of those L 2 functions on E which are Z-periodic. Hence 
this closed span cannot be equal to L 2 (E), so that the functions {1 / \J ' /i(E)e~ 27Tlkx \ 
k £ Z} do not form a frame for L 2 (E). Just as in Example 14.21 this implies that 
the set {D^Tktp \ m, k £ Z} cannot form a frame for L 2 (R) in this case. 

5. Multiresolution Analysis, Generalized Multiresolution Analysis, and 

the SNAG Theorem 

We now discuss multiresolution analyses and generalized multiresolution analy- 
ses. The concept of multiresolution analysis (commonly abbreviated by the acronym 
MRA) was originally developed by S. Mallat and Y. Meyer in the mid- 1980's 
( Ma, Me as a new way to construct interesting frames and wavelets. A priori, it 
would appear that their construction is not connected to the work of Stone and/or 
von Neumann. If one looks at the basic algorithm from a slightly different ap- 
proach, especially in the generalized multiresolution analysis case, it becomes clear 
that there is a connection. We first review the Stone-Naimark-Ambrose-Godement 
generalization of the Spectral Theorem for bounded self-adjoint operators due to 
D. Hilbert, E. Hellinger, and H. Halm. 

5.1. Unitary representations of locally compact abelian groups: the 
SNAG theorem. In 1930, building on earlier observations of H. Weyl, M. Stone 
noted the connection between unitary representations of R on a Hilbert space Tl, 
unbounded self-adjoint operators on H, and projection valued measures (resolutions 
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of the identity) on TL. To be more specific, he proved that associated to continuous 
unitary representation t \— > Ut of R on TL there is a unique (possibly unbounded) 
self-adjoint operator A defined on TL such that for every v, w G TL, Vt G R, 

< £/"*(«), W > = e 2mt<Av - w> . 

We express this in operator terms by writing 

U t = e 2witA , t G A. 

By the ordinary spectral Theorem for unbounded self-adjoint operators, we can 
write 

A = xdP r 



where P is a projection valued measure on (R, B), where B denotes the Borel subsets 
of R. Hence we can write 

U t = I e 2mtx dP x , Vt G R. 



By results relating projection valued measures to direct integral theory (see The- 
orem III. 13 of Fabec's book |Fap this implies that there exists a Borel measure v 
on (R, B) such that TL can be decomposed as a direct integral 



U = / [H x ]dv. 



The class of the measure v is unique. With respect to this direct integral decom- 
position of 7i, we can write 



(5.1) U t = / [e 27Tltx Id]dv(x), t G R, 

JR 

where Id denotes the identity operator on Ti x . Now let G denote the Pontryagin 
dual group of characters of a locally compact abelian group G. Note that because 
R = R via the correspondence x h- > e 2mx ' — \x-, Equation 15. II can be rewritten 
more generally as 

U t = / [ x (i) lA]dv{ x ), xei,teK, 

JR 

for a suitable Borel measure on R. 

This result was generalized to an arbitrary locally compact abelian groups by 
W. Ambrose [A|> M. Naimark |N] and R. Godement [SI in the early 1940's. Indeed 
G. Mackey calls the general result the "SNAG" theorem, for obvious acronymic 
reasons. 

Theorem 5.1. [SJ IT?l 1X1 1^] Let G be a locally compact abelian group, and 
let U be a continuous unitary representation of G on the Hilbert space TL. Let G 
denote the Pontryagin dual group of G. Then there is a Borel measure v on G such 
that Ti. can be decomposed as a direct integral 

H = / ffi [W x ]^( x ), 

JG 

and with respect to this decomposition we can write 

u 9 = C [x{g)Id\dv{ X ), X € G, g G G. 

JG 
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The measure class of v is unique. 

This result specialized to the case G = R gives Stone's Theorem [5]. 

In wavelet theory, this result will most often be applied to the case where 
G = Z", and Ti is a closed subspace of L 2 (R n ), usually denoted by Vo, which is 
invariant under the translation operators {T v : v £ Z n }. 

5.2. Multiresolution Analysis, and the construction of wavelet and 
frame bases. We give now the definition of Multiresolution Analysis due to S. 
Mallat and Y. Meyer. 

Definition 5.2. Let {Vi}i e % be a bisequence of closed subspace of L 2 (R"), 
and let A be an n x n dilation matrix, that is, a matrix all of whose eigenvalues 
have modulus greater than 1. We say that {Vijigz is multiresolution analysis 
(MRA) for dilation by A if: 

(i) ■ • • V-i C Vq Q V\ ■ ■ ■ (the Vi form a nested sequence of closed subspaces 
ofL 2 (R™)); 

(ii) D\{V Q ) = V h i e Z; 

(iii) U ieZ Vi = L 2 (R"), n igZ V5 = {0}; 

(iv) There exists 4> E L 2 (R), called a scaling function for dilation by A, such 
that 

{T v (<p) : u £ Z"} is an orthonormal basis for Vo- 

The last condition can be interpreted as meaning that Vq can be identified with 
Z 2 (Z n ), and hence by Pontryagin duality, with L 2 (T n ). Using this identification, we 
see that the representation of Z™ on Vo given by the operators {T v } is equivalent to 
the regular representation of Z™ on l 2 (Z n ). From the direct integral point of view, 
using the SNAG theorem just discussed, this is equivalent to saying that the scaling 
function gives us a spectral decomposition of the translation operators T v , v € Z™, 
acting on of the space Vo- Let us set n = 1 to illuminate this idea more clearly. The 
measure on Z = T is guaranteed by the SNAG theorem turns out to be the Haar 
measure on T. Define a map J : Vq — > L 2 (T) by 

J(]Tc fc T fc (0))(x) = ]Tc fce 2 ^. 

Then using elementary facts about orthonormal bases in Hilbert spaces, it is easy 
to see J is a unitary isomorphism, and that for every f £ Vo, 

J(T(/))(x) = e 2 ™J(f)(x). 

We will return to this construction when discussing the generalized multiresolution 
analyses of L. Baggett, H. Medina, and K. Merrill BMM . 

We first view MRA's in the frequency domain by means of the Fourier trans- 
form. 

Definition 5.3. Let {V^} ie z be a bisequence of closed subspace of L 2 (R"). We 
say that {Vi]i £ z is a multiresolution analysis in the frequency domain for 
dilation by A if: 

(i) ■ • • V-i C Vq C Vi ■ ■ ■ (the Vi form a nested sequence of closed subspaces 
of J L 2 (R")); 

(ii) (DlYi Vo) = V u i G Z; 

(iii) U teZ V = i 2 (R"), n ie zV = {0}; 
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(iv) There exists $ € L 2 (R n ), a scaling function in frequency domain, such 
that 

{T V (Q) : v £ Z™}, is an orthonormal basis for Vq. 

Of course one moves between MRA's and MRA's in the frequency domain by 
setting 

Vi = F(Vi), V = F-\Vi) 
Example 5.4. Let n — 1 and A = (2). Then 



% - 1 , if x ^ 0, 
1 if x = 



is a scaling function in the frequency domain for dilation by 2. 

Let A be a n x n integral dilation matrix. Given a MRA for dilation by A, 
there is an algorithm using the scaling function $ for constructing a wavelet family 
in the frequency domain for dilation by N. In this algorithm, the scaling function 
$ is extremely important, since from <I> one constructs Vq using the translation 
operators T v , and then the other subspaces Vi using the dilation operator D, and 
finally the wavelet family from looking at Wo = V\ Q Vq. 

In order to construct it helps to have a low-pass filter m for dilation by 
A. For simplicity of exposition we let n = 1 and let the dilation matrix be N where 
N is a positive integer greater than 1. 

Definition 5.5. Let N be a positive integer greater than 1. A low-pass filter 
mo for dilation by TV is a Z-periodic function mo : R™ — > C which satisfies the 
following conditions: 

(i) too(O) = v^; 

(ii) E;=o 1 K(z + Z/AO| 2 = iVa.e. 

(iii) Too is Holder continuous at and is non-zero in a sufficiently large neigh- 
borhood of ("Cohen's condition", c.f. [HolID^ 

Given a low-pass filter, one can use it to construct the scaling function in the 
frequency domain as follows: 

Theorem 5.6. [Mai IMel lD2] Given a low -pass filter too for dilation by N 
which satisfies Cohen's condition, then 

converges a.e. and in L 2 (M.) to a scaling function in the frequency domain for 
dilation by N, which can be used to construct an orthonormal wavelet family for 
dilation by N. 

We briefly discuss how one uses the scaling function, the low-pass filter, and 
other filters called high-pass filters, to construct the wavelet family. 

Definition 5.7. Let n—1 and let the dilation matrix A now be the positive 
integer N, where N > 1, and let too be a low-pass filter for dilation by N. A set of 
essentially bounded measurable Z periodic functions mi, 7712, ■ ■ ■ , m^-i are called 
high-pass filters associated to too, if the N x N matrix 

C i{X+ J /N) )o<i,j<N-i G U(N,L°°(R/Z)). 
v N 
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The following result is due to Y. Meyer and S. Mallat, and in the normalized 
tight frame case to W. Lawton. 

Theorem 5.8. |MaL ICTel HI] Let N be a positive integer greater than 1, let m^ 
be a low-pass filter for dilation by N, and let $ be the scaling function constructed 
from tuq above. If mi, m%, ■ ■ ■ , mjsr-i are high-pass filters associated to mo, then 

{*i = 5(mi$), W 2 = D{m 2 <S>), • • • , *jy-i = 5(mAr_i$)} 

is a orthonormal wavelet family in the frequency domain for dilation by N. If Co- 
hen's condition is not satisfied, the {^k} '$ still form a NTF wavelet family. 

The above discussion indicates how to obtain MRA's from filter functions. It 
may not be immediately apparent how filter functions are related to the work of 
von Neumann, but a very nice paper by L. Baggett, A. Carey, W. Moran, and P. 
Orring BCMQ] gives an exposition of wavelet theory from the point of view of 
von Neumann algebras, and shows how one can obtain filter functions from MRA's 
by using a result on cancellation for finite von Neumann algebras. The result 
of Baggett, Carey, Moran and Orring employed von Neumann algebras to study 
certain lattices F of locally compact abelian groups G having a faithful unitary 
representation U on a Hilbert space TL via translation. They also assumed the 
existence of a generalized "dilation operator" 8 acting on Ti with the property that 
{S~ 1 U 1 8 : 7 € r} was a subgroup of finite index in {U 1 : 7 € T}. Here 
for simplicity r = Z, G = M, H = L 2 (R), and Z acts on L 2 (R) by translation 
operators. We also assume that the dilation operator 6 is the standard dilation Djy 
corresponding to the positive integer N > 1, which we denote here by D. We refer 
the reader to |BCMO] for the general set-up; we only note that the main theme 
of the argument carries over to the general case. 

Let {Vj}j£% C L 2 (R) be an ordinary multiresolution analysis for dilation by 
N. Then Vq is invariant under translations {T v : v € Z}, and by definition, the 
corresponding representation of Z on Vq is equivalent to the regular representation 
of Z on Z 2 (Z). Let W*(Z) denote the commutative von Neumann subalgebra of 
B(L 2 (R)) generated by {T v : v S Z}. Then V x = D(V Q ) is invariant under the 
action of W*(Z) with invariant subspace Vq, and the multiresolution theory tells 
us that 

Vx = Va® Wo, 

where Wq is the wavelet space, which is an invariant subspace for the von Neumann 
algebra W*(Z). Using cancellation properties for finite von Neumann algebras (c.f. 
|Dx| . Section 3.2.3 Proposition 6), they are able to show the following result: 

Theorem 5.9. BCMO The wavelet space Wq can be decomposed into N — 
1 invariant orthogonal subspaces {W^} 1 ^ 1 , each invariant under the action of 
W*(7j). For each i, the action ofl on the representation o/Z on Wq corresponding 
to the action ofW*(7i) is equivalent to the regular representation of Z on Z 2 (Z). 

By definition Vq contains an element cj> such that {T v ((f>)} ve i forms an orthonor- 
mal basis for Vq. Viewing this in the Fourier domain, 



{T v (0)} ve z = {T„($)} uez 
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forms an orthonormal basis for Vq = T(Vq), where $ = <j>. Since $ G Vi 
D(Vo), we can write 

*(s) - D[F(J2ao,vT v (<t>))](x) = a ,„ e 2irira )*](a;), 

where the infinite sums converge in the relevant Hilbert space norm. Writing 

m (x) = Y, a ^e 2mvx , 

we get 

= D(m $)(x) = -Lm (^)$(^), 



which is exactly the refinement equation in the frequency domain. The function 
mo turns out to be the low-pass filter in this situation. 

Similarly, each subspace Wq, 1 < i < N — 1, contains an element ipi such that 
{T v (ipi)} ve i forms an orthonormal basis for Wq. Again we have 

{T v {$i)}ve-z = {T v (^t)}vez 

forming an orthonormal basis for Wq = T{Wq), where ^ — ipi. Since \&j € V\ = 
D(V ), i G {1, • • -N — 1}, we see that & G Vq with ipi = D(£i). Then as before, 



Writing 



we obtain 



mi (x) = Va v e 2 ™ OT , 1 < N-l, 



^i(x) = D(mi$)(x), 
and the family {m 1; to.2, ••• , tojv_i} are the high-pass filter functions associated to 
toq- Standard arguments from wavelet theory (c.f. |Str| ) show that the filter func- 
tions {too, mi, to-2, • • • ,mN~i} constructed in this fashion satisfy the standard 
high-pass orthogonality conditions as condensed in Definitions 15 . 51 and 15 . 71 

An approach to finding continuous high-pass filter functions corresponding to a 
given low-pass filter function by using a C* -module point of view, which has features 
in common with the approach in |BCMO] , can be found in a recent article by the 
author and M. Rieffel |PR) . 

5.3. Generalized multiresolution analyses and relationships to the 
SNAG Theorem, and current research. In the frequency domain version of 
MRA for n = 1, Condition (iv) implies that the representation of Z on the initial 
space Vq is equivalent to the regular representation of Z on l 2 (1). This means 
that when we use the Spectral Theorem to rewrite this representation as a direct 
integral on Vq, there is a measure rfv on Z = T = R/Z, which we parametrize 
by x G [0, 1), such that 

r® 

V = / [C] dv(x) 

and T v = M e 2nivx dv(x). It is clear from our Fourier transform of T that dv{x) 
is just Haar measure dx on T. The use of the SNAG Theorem comes when the 
wavelet theory is extended to translations T v , v G Z", acting on L 2 (R"), and 
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more general dilations Da corresponding tonxn integer dilation matrices A. In 
the ordinary MR A case, we get a somewhat analogous situation to the case when 
n = 1. Things become more complicated when condition (iv) of the definition of an 
MRA is not satisfied. 

We give an example of a wavelet which does not come from an MRA. 

Example 5.10. Consider the wavelet in the frequency domain coming from 
the wavelet set 

r 16 „ s r 1 2 X r 2 1. . 16, 
E = [-y,-2)U[--,--)U[-,-)U[2, y ). 

Then xe — j is known as the Journe wavelet, and is a wavelet for dilation by 
2 in the frequency domain that does not come from an MRA . However, one can 
still form the nested sequences of Hilbert spaces 

V = sp^{D j T v {T- l {^j)) | v G Z, j < i] 

and it is possible to show: 

(i) ■ • • V-i C Vq C Vi • • • (the V{ form a nested sequence of closed subspaces 
of£ 2 (R)); 

(ii) D^Vp) = V u i G Z; 

(iii) U iEZ Vi = i 2 (M), n ieZ V t = {0}; 

(iv) Vo is invariant under all powers of T. 

The situation when one has a wavelet family and NTF wavelet frame which 
does not come from an MRA has been approached in various ways in independent 
works of L. Baggett, H. Medina and K. Merrill jBMMj. J. Benedetto and S. Li 
|BL| . and D. Han, D Larson, M. Papadakis, and T. Stavropoulos HLPS . Wc 
discuss here the approach due Baggett, Medina and Merrill. 

Definition 5.11. Let A be a n x n integral dilation matrix. A generalized 
multiresolution analysis (GMRA) for dilation by A is a sequence of closed sub- 
spaces {Vi}i£% of L 2 (M. n ) which satisfy the following conditions: 

(i) Vi C V i+ i for all i G Z; 

(ii) D A {Vi) = Vi+ U VzeZ; 

(iii) U ieZ Vi - i 2 (M"), r\ezVi = {0}; 

(iv) Vq is invariant under all of the operators {T v : o £ Z™}. 

One can reformulate the above conditions in the frequency domain if desired, 
just as in Definition 15. 31 

Baggett, Merrill and Medina showed given a wavelet family, one can correspond 
to it a GMRA, as in Example l5.10l Given a GMRA not associated to a wavelet, they 
also were able to develop several characteristic invariants associated to a GMRA 
using spectral theory. They generalized the unitary operator J defined from Vo to 
L 2 (T) in the classical MRA case to the GMRA case in the following theorem. 

Theorem 5.12. BMM Given a GMRA in L 2 (K n ) corresponding to dilation 
Da by an integral dilation matrix A and the transform of integer translation T v , 
there is a unique sequence of Borel subsets Si 2 S2 2 • • • of T™ and a unitary 
operator J : Vq — > @j L 2 (Sj) such that 

[J{T v {f)Ux) = e 2m [J(fMx), /eL 2 (T"), v G Z». 
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We briefly discuss the construction of J and the sets {Sj} from the point 
of view of the SNAG Theorem. The representation T of Z™ on Vq can be de- 
composed as a direct integral. With respect to this direct integral we can write 
Vo = }^ n [Hx\dv(x), where v is a Borcl measure on Z™ = T", and then by the 
SNAG theorem, with respect to this decomposition we can write 

T v = e 2m ld]dis(x), v e Z n . 

Baggett Merrill and Medina first proved that the class of the measure v coming 
from the SNAG theorem is exactly the class of Haar measure on T™, and then set 

Sj = {x € TP dim(?4) > j}. 

Baggett, Merrill and Medina called the function = J2jXs'j( x ) defined 

on T the multiplicity function corresponding to the GMRA {T^}i S z, because 
n(x) indicates the multiplicity of the character e 27rzx ' e Z" = T™ in the above 
decomposition of the representation T n , i.e., 

fi(x) = dim(H x ), x G T". 

For simplicity, we specialize to the case where n = 1 and A is dilation by a 
positive integer N > 1. Baggett, Merrill and Medina showed that /j, satisfies 

JV-i x + l 

mO) < ^("at) a - e - 

1=0 

If /i is essentially bounded, with c = ess sup/x(x), it is possible to define the conju- 
gate multiplicity function 

1=0 

By definition, /i and // satisfy the so-called "consistency equation" : 

N-i , 

(5.2) (j,(x) + = Y M(-^-)- 

Baggett, J. Courter and Merrill in BCM then generalized the Mallat and 
Meyer algorithm for constructing wavelets from filters to this GMRA setting. They 
first generalized the concept of low and high pass filters. Given an integer valued 
function \x (essentially bounded by c on T and satisfying technical conditions out- 
lined in BM they defined "generalized conjugate mirror filters" , which correspond 
to low-pass filters in the classical case, to be functions {hij}i< ij < c , where each hij 
is supported on Sj. Similarly, they defined "complementary conjugate mirror fil- 
ters," an analogue of classical high pass filters, to be functions {gk,j}i< k <d, i< j <o 
where each g^j is supported on Sj. 

The functions {hij} and {gk.j} satisfy orthogonality conditions which are mod- 
ified versions of the orthogonality conditions satisfied by classical low and high-pass 
filters; see |BCM| for precise details. They built examples of functions gtj and dis- 
satisfying these conditions by using an explicit algorithm, and then modified the 
results obtained to get examples of generalized filters, mimicking the way examples 
of filters are obtained in the classical case BCM . 
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Under appropriate conditions on these generalized filter functions, Baggett, 
Courter and Merrill then used them to construct a finite tight frame wavelet family 

{*!,-•• CI 2 (I). _ 

They first built generalized scaling functions {t&i, <I?2, • • • , ^c} Q Vo using an 
infinite product construction involving dilates of a matrix with periodizations of 
{hi,} as entries (c.f. Theorem 3.4 BCM). The appear as the entries in 

the first column in the infinite product matrix. 

Given the above notation and construction, they then have: 

Theorem 5.13. |BCM| Let {(ft»,j}i< i,j <c and {gk,j}i< fe <d, i< j <c be gen- 
eralized filter functions associated to the multiplicity function /i, that satisfy appro- 
priate conditions and let {$i,$2,-- - , $ c } C Vq be generalized scaling functions 
constructed as described above. Setting 

c 

* fc = D(J2 9k,j$j), 1 < k <d, 

3=1 

the {^k | 1 < k < d} form a NTF wavelet family for dilation by N. 

Current research of the author, Baggett, Merrill, and P. Jorgensen involves a 
generalization of certain results of Baggett, Courter and Merrill. Bratteli and Jor- 
gensen have found a group (called the loop group) which acts freely and transitively 
on Lipi filter systems {mo, mi, • • • , TOjv_i} for dilation by N (such a collection of 
filter functions is called an m-system by Bratteli and Jorgensen) , and have shown 
that each Lip\ m-system can be used to construct a NTF wavelet family with N — 1 
elements. Just as m systems have been shown to give rise to functions on T with 
values in N x N unitary matrices, we have shown that each generalized filter system 
(hij) and (gk,j) coming from a GMRA of Baggett, Courter and Merrill gives rise to 
unitary-matrix valued functions on T, whose dimensions now vary with We 
have also found an analogue of the loop-group which acts on the generalized filter 
functions {hij) and {gu.j) to give new functions satisfying the same orthogonal- 
ity relations |BJMP| . We conjecture that operator algebraic methods in addition 
to the results of BJMP to will allow us to generalize Theorem 5.13 of Baggett, 
Courter and Merrill stated above. 

As we have seen, operator algebras and direct integrals have played key roles 
in abstract wavelet theory, and we hope they will help in our new research project 
as well. 
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